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ABSTRACT 

The  theorems  of  limit  analysis  are  not  valid  if  part  of 
the  boundary  is  pressed- against  a rigid  surface  and  constrained 
by  coulomb  friction.  In  reference  to  a problem  of  inverted 
extrusion,  we  show  how  conventional  limit  analysis  can  be 
extended  to  obtain  upper  and  lower  bounds  on  the  force  necessary 
for  extrusion. 
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LIMIT  ANALYSIS  AND  COULOMB  FRICTION 


When  part  of  a rigid/perfectly-plastic  body  is  pressed 
against  a surface  and  constrained  by  coulomb  friction,  then 
the  normal  velocity  component  will  be  zero.  The  tangential 
component  will  also  be  zero  unless  the  shear  traction  stress 
is  equal  to  uN  where  y is  the  coefficient  of  friction  and  N 
is  the  compressive  normal  traction  or  is  equal  to  the  shear 
yield  stress  k.  Thus 

0 £ T <_  min  (k,  yN) 

> 0 (1) 

IF  > 0 then  T = min(k,  yN) 

where  T and  are  positive  in  opposing  directions.  Equations  (1) 
do  not  explicity  prescribe  either  stress  or  velocity,  and  hence 
the  conventional  theorems  of  limit  analysis  are  not  generally  valid. 


As  an  excunple  of  a case  where  this  phenomenon  occurs,  we 
consider  inverted  extrusion  through  a tapered  die  with  coulomb 
friction.  Fig.  1.  The  yield-point  load  F is  that  value  of  the 
applied  force  which  will  just  cause  the  billet  to  move  through 
the  die  at  an  arbitrarily  slow  speed.  We  consider  the  case  of 
plane  strain;  a similar  result  obviously  applies  to  rotational 
symmetry . 

We  shall  assume  that  at  yield-point  motion  takes  place  all 
along  the  die  face  AB.  Then  the  kinematic  and  static  boundary 
conditions  are,  respectively 


On  ACC "A':  V = V, , V = 0 
X i.  y 

On  AB  and  A'B':  = 0,  > 0 


Viti/t2, 


(2) 


On  BDD'B"; 


0 


2 


and 


On  AB  and  A'B': 


T = (1/2) (a  -a  ) sin  2a  - x cos  2a  (3a) 

y X xy 

N = - (1/2) (a  +a  )-(l/2) (a  -o  )cos  2a  - x sin  2a  (3b) 
y X y X xy 


T = min  (k,  yN)  N > 0 


(3c) 


Over-all  horizontal  equilibrium  of  the  billet  shows  that 

S, 


tj^P  = / (T  cos  a + N sin  a)ds 


(4) 


Since  it  follows  from  Eq.  (3c)  that  N ^ T/u,  Eq.  (4)  leads  to 
the  inequality 

^1 

/ Tds  ^ tj^F(cos  a + sin  a/u)~^ 


(5) 


We  define  a statically  admissible  stress  field  as  stresses 
which  satisfy  internal  equilibrium,  the  yield  condition  and 
boundary  conditions  (3) , together  with  a load  F®  defined  by 
(4) . The  principle  of  virtual  work  then  shows  that 


1 / a.®  e.  . =*  t,V,F®  - f T®V.ds 

2 - 1]  1]  11  t 


(6a) 


1 ^ 

T / “ t.V  F - f ^ TV  ds  (6b) 

z - ij  IJ  X j.  u 

" S 

^2 

Subtracting  (6a)  from  (6b)  we  obtain 


tj^Vj^(F-F®)  - f ^(T-T®)V^ds  = j ^ij  1 0 C7) 


where  the  last  step  follows  from  Drucker's  postulate  [1]  for 
a perfectly-plastic  material. 


V 
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Since  T® , T,  and  are  each  non-negative  at  each  point, 
Eq.  (7)  implies  the  weaker  restriction 

Si 

t^V^F  - 0 1 t^V^F»  - / TMs  (8) 


Finally,  we  make  the  physically  reasonable  assximption  that  since 
the  frictional  force  is  restraining  the  material  at  the  boundary,  the 
horizontal  component  of  the  boundary  velocity  will  not  be  greater 
than  the  average  horizontal  component  across  a section  of  the 
billet.  Thus 

t max  X max  — i 1 z 

so  that  we  obtain  the  lower  bound 

^1 

F ^ F®  - / T®  ds/(tj  COS  a)  (10) 

For  an  upper  bound  we  start  with  any  incompressible 
velocity  field  V^*  which  satisfies  (3) , define  a strain  field 
by  e. .*  = (1/2) (V.  .*  + V.  .*)  and  a stress  field  o . by  the 

normality  flow  rule,  and  define  the  dissipation 

D*  = i ;o.  *e. .*  (11) 

Z 1 j i J 

The  principle  of  virtual  work  shows  that 
1 • ^1 


Subtracting  both  sides  of  (12)  from  D*  and  using  Drucker's 
postulate  [1]  we  obtain 


D*  - + / TV^*ds  = j i ° 


whence 


tiViF  - (V^*), 


/ Tds  < D* 


Finally,  we  use  (5)  to  obtain  the  upper  bound 


(V^*)  /V, 


< ^ r _ max 

— ^ cos  a + sin  a 


Alternative  bounds  can  be  obtained  by  considering  the 
extreme  cases  of  perfectly  smooth  and  perfectly  rough  dies. 
Obviously  the  yield-point  load  for  the  coulomb  die  will  lie 
between  these  two  extremes.  Therefore,  any  upper  bound  on  the 
rough  die  or  any  lower  bound  on  the  smooth  one  will  also  be  a 
corresponding  bound  on  the  coulomb  die. 

A perfectly  rough  die  can  permit  motion  only  when  T = k. 
Starting  with  any  kinematically  admissible  velocity  field 
we  proceed  as  above  but  replace  T by  k in  (12) , thus  leading 
to  the  alternative  upper  bound 

®1 

F < (1/tj^Vj^)  [D*  - k f ^ V^*ds]  (16) 

For  the  prefectly  smooth  punch,  no  shearing  force  can  be 
supported  along  AB,  so  that  Eq.  (3c)  is  replaced  by 

On  AB  and  A'B':  T » 0 (17) 

Denoting  a statically  admissible  field  under  this  changed 
boundary  condition  by  we  obtain  the  lower  bound  of 

conventional  limit  analysis 

F > F®®  (18) 

In  closing,  we  emphasize  two  restrictions  on  the  results 
herein.  In  conventional  limit  analysis,  the  minimum  upper 
bound  is  equal  to  the  maximum  lower  bound  and  this  common  value 
is  the  exact  yield-point  load.  Here,  however,  the  maximum 


mM 


of  all  (10)  or  (18)  will  generally  be  distinctly  less  than  the 
minimum  of  all  (15)  or  (16) , so  that  we  can  never  do  better 


than  define  a range  within  which  the  yield-point  load  may  lie. 
Secondly,  our  results  depend  upon  the  assumption  that  is 
strictly  positive  and  bounded  by  (9) . Although  this  assumption 
is  physically  reasonable,  we  have  not  proved  its  necessity. 
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